Cavity photon counting: ab-initio derivation of the quantum 
jump superoperators and comparison of the existing models 

Teppo HayrynenI*] Jani Oksanen, and Jukka Tulkki 
Department of Biomedical Engineering and Computational Science, 
Helsinki University of Technology, P.O. Box 9203, FIN-02015 HUT, Finland 

(Dated: August 12, 2008) 

Abstract 

Time development of electromagnetic fields in closed cavities under continuous detection of 
photons continues to be a subject of confusing controversy. Recently Dodonov et al. [Phys. Rev. 
^ ' A, 75, 013806, 2007] argued that their model of quantum superoperators (E model) invalidates 
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some of the predictions of the previously introduced photon counting model of Srinivas and Davies 
[J. Mod. Optic. 28, 981, 1981] (SD model). Both the SD and the E models are based on 



• two postulated quantum jump superoperators: (1) the one-count operator corresponding to the 



absorption of a single photon and (2) the no-count operator. In this work we develop a stochastic 
^ i difference equation that describes the dissipative coupling of the cavity field and the detector. The 

difference equation is based on non-perturbative treatment of the cavity-detector coupling. In spite 
of being non-integrable due to the coupling of the detector with an external reservoir it can be 
, used to derive the exact forms of the quantum jump superoperators. When applied to a particular 

' photon counting measurement our theory gives predictions identical with those of the SD model 

OO ' which should be considered a non-perturbative and ab-initio result. It is pointed out that available 

o . 

^0 . experimental results coincide with the results given by the ab-initio SD model. We summarize 

some of the key characteristics of cavity fields and photon counting processes to demonstrate that 
^ ' the results given by the SD model are consistent with the principles of quantum mechanics while 

those given by the E model are not. 
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I. INTRODUCTION 



The theory of cavity photon counting and in particular the expected evolution of the 
cavity field after detection of a sequence of photons continues to be a subject of confusing 
controversy as shown in a recent report by Dodonov et al. [1]. So far no experimental data on 
the evolution of the cavity fields are available to resolve the controversy . The experiments 
themselves should be feasible at present though. This has prompted us to study the relation 
between the photon counting models. In this work we consider a system formed by the 
cavity field, the detector, and a (infinite) reservoir that is coupled to the detector. 

The first quantum mechanical models of photon counting were formulated for an exper- 
iment where a light beam entered the detector and the unabsorbed photons were able to 

n n 

escape p^3|. Those models were based on the quantum version of the classical Mandel's 
formula 4|- In these models there was thus no need to consider the measurement back 
action. However, in the case of a cavity field the time integrated effect of the measurement 
on the electromagnetic field will be large and therefore the measurement back action cannot 
be ignored. We will show that the analysis of photon detection measurements can be carried 
out starting from the interaction Hamiltonian of the field and the detector. It is pointed 
out that for the photon counting analysis there is no need of introducing new postulates 
beyond the well established theoretical foundations of cavity quantum electrodynamics. It 
;urns out that our non-perturbative approach gives results which agree with the SD model 
4] while the E model jll, [sl, [g, Q] is shown to be inconsistent. Furthermore, we show that 
the measured data of the second order coherence degrees of thermal and coherent fields 81] 
agree with the SD model. 



II. PHOTON COUNTING BASED ON OPEN SYSTEM MASTER EQUATION 

The SD and E photon counting models were originally introduced using phenomenolog- 
ical arguments or postulates. The first rigorous approach to cavity photon counting was 
formulated by Imato et al. 9|. They used a homogeneous atomic beam as a tool for includ- 
ing the dissipative terms in the Jaynes-Cummings Hamiltonian. The dissipation of the field 
was modeled using a sequence of infinitesimal perturbations caused by the atoms. 

In this work we show that the atomic beam model is not needed for derivation of the 
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theory of cavity photon counting. Instead we assume that the whole system consists of three 
parts: the cavity, the detector and a (infinite) reservoir. The details of the coupling between 
detector and reservoir arc unknown but it is assumed that this coupling is strong enough 
so that the detector is most of the time in its ground state \g). The time the detector will 
spend in its lowest excited state |e) is so short that the probability of a cavity photon to 
interact with the detector in the excited state is vanishingly small. Accordingly in any short 
time interval only absorption of one photon or no photons are possible i.e. the one-count or 
no-count events, respectively. The photoabsorption processes are Markovian processes i.e. 
the system is memoryless and the absorption probability depends only on the current state 
of the system, not on the past states. Furthermore, we assume that the detector doesn't 
return to the ground state by emitting the photon back to the cavity. 

The optical cavity is assumed ideal and thus the only dissipative mechanism included is 
the absorption of photons by the detector. Furthermore, because the detector is assumed 
to return back to the ground state immediately after the absorption of a photon we do 
not consider the detector saturation effects. If we assume that the photons escape from 
the cavity to a detector which is, for example, a photomultiplier tube, the photons would 
dissipate from the cavity to the photomultiplier tube even if the tube was be momentarily 
switched off. Thus we can neglect the dead time of the detector. 

We assume that the interaction of the field and the detector is described by the Jaynes- 
Cummings Hamiltonian in the rotating wave approximation 

Hi^hn(a^\g){e\ + a\e){g\y (1) 

Furthermore, the field Hamiltonian is hud^d, the detector Hamiltonian is hiu\e){e\, and the 
initial detector-field density operator is 

oo 

Pii^)= ^ Pn,n'\n,g){n',g\. (2) 

n,n'=0 

Wc next find an expression for the differential change of the detector-field density operator 
dpi — pi{dt) —pi{0). At the starting point of the interval [0, dt] the detector is in the ground 
state and the probabihty of the detector being excited to |e) is small. Therefore during [0, dt] 
the coupling of the detector with the reservoir can be neglected since the reservoir only causes 
relaxation of the detector from the excited state to the ground state. Thus the field detector 
subsystem can be considered closed within [0,dt]. The differential change of the density 
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operator can then be obtained non-perturbatively using the Taylor series 



piidt) - p,(0) 



dpiit) 



dt 



dt + 



1 d'piit) 



2 dt^ 



(dtr + .... (3) 



i=0 

As shown in detail in the Appendix I A 11 the Taylor expansion gives (see Appendix lA 1 [ ) 

P'^'^^^-^P'^^^ = -^[Hj,pm] + ^ {2HM0)Hi - {HjHiJim) . (4) 

Note that this result is exact in the limit of infinitesimal dt. From equations ([T]) 
and we see that Hipj{0) = JT'^J2nn'Pn,n'Vn\n — l,e){n',g\ = hQa\e) {g\pj{0), that 
HiHiPiiO) = {hn)^J2n,n'Pn,n'n\n,g){n',g\ = {hnya^api{0), and that Hipi{0)Hi = 
(hiiy J2n n' Pn,n'Vnn'\n — l,e){n' — l,e| = {hil)'^a\e) {g\pi{0)\g) {e\a^ . Using these relations 
and replacing the initial time t = with t, equation (jlj) can be written as 

dp lit) ^ Piit + dt) - Piit) 
dt dt 

= -iVL^^pn,n'it){^/n\n - l,e){n',g\ - \/n'\n,g){n' - 1, e| j 

71, n' 

' Ta\e){g\pi(t)\g){e\a^ -{a^arpiit)}]. (5) 



Equation ([5]) can not be solved by simple integration since it has been assumed that in 
the beginning of each time interval [t, t + dt] the detector is at ground state and therefore 
the coupling between the the detector and the reservoir has been neglected. With these 
assumptions equation ([5]) describes the field-detector system as an open quantum system 
dissipating its energy to an infinite reservoir represented by the detector returning to the 
ground state infinitely fast. In the following we are interested in the evolution of the field 
density operator. This is obtained from equation ^ by calculating the trace over the 
detector states pf = {g\pi\g) + (e|p/|e). Setting dpj = pj(t + dt) — p/(t) and assuming that 
the detector is at ground state at t, the reduced density operator i.e. the density operator 
of the field can be written as 

dp/(t) _ ^'^tf2apf(^t)a^ -{a^ajfit)}], (6) 



dt 2 

Note that in contrast to equation ([5]), equation ([6]) may be solved by integration for 
our detector model, since it no longer includes the detector states. Substituting p/(t) = 
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Yl'^,n'=oPn,n'(t)\n){n'\ iiito equation © gives 

dp/W dpn,n'{t) , 

dt ^ dt ' ' 

n,n'=0 



(oo CXD , V 
^ ^ ^ _l ly^ \ 
2J yrm'Pn,„/(t)|ra - - 1| - 2J ^ Pn,n'(^)l^)(^'l ) • 



(7) 

' n,n'=0 n,ri'=0 

In order to relate the constant f2^dt to the couphng strength of the field-detector we solve the 
diff'erential equation for the probabilities p„,„ i.e. the master equation of the photon numbers. 
Taking the diagonal matrix elements (n| ■ |n) of equation ([7]) and denoting p„ = gives 

dpn(t) _ d{n\pj{t)\n) 
dt ~ dt 

= nMt((n + - (8) 

Thus, the detector only absorbs photons and do not emit them to the cavity. The constant 
A = n^dt is the probability per unit time that a photon will be absorbed by the detector 
i.e. the rate of absorption. 

The physical interpretation of ^dt becomes particularly transparent if one considers the 
absorption of a single photon by the detector. Let the amplitude of the detector ground 
state + one photon field state (|(yf, 1)) be given by Cgit) and the amplitude of the excited 
state (|e, 0)) by Ce(t) with initial conditions Cc,(0) = 1 and Ce(0) = 0, respectively. The 
differential change of Ce(t) can be then obtained from the time dependent Schrodinger 
equation. It is given by CJydf) = —iVe^gdt/h, where V^^g = {e,0\Hj\g, 1) = hQ is the dipole 
amplitude. The probability that the photon has been absorbed is |Ce(dt)p = | Vg^^^dt^/^^. 
Thus the absorption probability per unit time is given by \Ce{dt)\'^/dt = \Ve,g\'^dt/ti? = Q'^dt. 
This is in accordance with equation (JHl). 

III. ONE-COUNT AND NO-COUNT SUPEROPERATORS DERIVED FROM 
MASTER EQUATION 

The density operator of the system after infinitesimal time dt is obtained using equation 
(IS]) and relation pi{t + dt) = pi{t) + ^^^dt which gives 



pj{t + dt) = pi{t) + 



iVL'^Pn,n'it) [Vn\n- l,e){n',g\ - Vn'\n,g){n' - l,e|^ 

dt. (9) 



+^—{2a\e){g\pi{t)\g){e\a^ - {a^a,pi{t)} 



After measuring a photoabsorption process the detector is in the excited state |e) so the part 
of the field density operator corresponding to the one-count event is given by {e\pi(t + dt)\e). 
If a photon is not absorbed the detector stays in the ground state and the part of the field 
density operator corresponding to the no-count-event is given by {g\pi(t + dt)\g). We obtain 

{e\pi{t + dt)\e) = \apf{t)a^dt (10) 
{g\pi{t + dt)\g) = pf{t) - ^ {a^apf{t)+pf{t)a^a) dt. (11) 

Therefore, after the sudden absorption of a single photon the one-count operator is simply 
given by equation (fTOl) 

Jpfit) = \apf{t)a). (12) 

In contrast to the abrupt one-count process the no-count process changes the field density 
operator in a continuous manner. Therefore it is necessary to calculate the no-count operator 
in a finite time interval [t, t + r] between two absorption events. The density operator is 
changed by the no-count event by the amount dpf{t) = —^{a)a,pf(t)}dt, obtained from 
equation (fTTl) . whose solution in time interval [t, t-|-r] (here r is not necessarily infinitesimal) 
gives the no-count operator 

SrPfit) = e-^'''-pf{t)e~^'''\ (13) 

The probabilities of absorbing a photon and not absorbing a photon are obtained 
from equations (|T2l) and (fT3ll . respectively, by calculating Trace{ Jp/}/Trace{p/} and 
Trace{S'T-p/}/Trace{p/}. Equation (fT2l) gives the probability of the abrupt absorption pro- 
cess at [t,t + dt] 

Pone-count{t, t + dt) = \n{t)dt. (14) 

The no-count probability at time interval [t,t + r] is obtained calculating the trace of (fT3ll . 
Using the series expansion of the exponential terms gives (see Appendix IA4I) 

oo 

Pno-countit, t + t) = ^ e-^''^Pn{t). (15) 
n=0 

Note that for infinitesimal time interval t = dt one obtains Yl'^=o^~^^'^'^Pn{t) = Yl'^=oi^ " 
\ndt)pn(t) = 1 — \n(t)dt i.e. the sum of the one-count and no-count probabilities is unity. 
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The change in the expectation value of the number of photons in a small time interval dt 
is obtained from equation © 

d|t) ^ f2n^ = xf:ni-np4t) + in + l)p^^,it)) 

n=0 n=0 
oo 

= -X^npnit) = -Xn{t) (16) 

n=0 

which gives 

n{t) = n{0)e-^'. (17) 

Thus the expectation value of the number of photons decays exponentially in time. 

Note the difference between our model and the model of Imato et al [9]. By using 
the perturbation theory, they calculated how an atom perturbs the cavity field during an 
infinitesimal interaction time and thereby dissipates the photons from the cavity. In contrast 
we have derived, starting from the time- dependent Schrodinger equation, a master equation 
for a system where the cavity field is coupled with a detector which relaxate to a infinite 
reservoir. However, both models leads to the SD model when the detector states are reduced. 



IV. COMPARISON OF PHOTON COUNTING MODELS 

The SD and the E models share four main postulates: 

(1) The absorption of photons the detector takes place as instantaneous events represented 
by the one-count superoperator Xipj(t) = 7^Apj(i(:)A^ (7^ is defined in postulate 4). The 
one-count operator is a quantum jump superoperator i.e. the density operator jumps from 
p/(t) to 7A^P/(t)^^ in infinitesimal time interval \t,t + dt]. In the SD model A = a. The 
E model is obtained from the SD model by replacing the well known bosonic annihilation a 
and creation a) operators by the normalized operators A = [a)a + l)~^/^a (denoted by E 
below) and A'^ = d^(a^d + 1)~^/^ (denoted by E"^ below), respectively |5|]. These normalized 
operators obey the relations E\0) = 0, E\n > 0) = \n — 1) and E^n) = |?t, + 1). The 
probability of the one-count at [t,t + dt] is 7ATrace{A/)/y4^}dt. 

(2) Between the counts the density operator evolves according to the no-count superoperator 
Srpf = e^^^^^pf{t)e^^^^/^, where Ya =-iHo - ^h-fAA^A and Hq = hua^d. Here r is not 
necessarily differential. 

(3) After measuring an event corresponding to the operator O {O = J ot S), the density 
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operator is Pf(t^) = Opf(t)/TTa.ce{Opf(t)}. 

(4) Furthermore, the couphng between the detector and the field is parameterized using a 
model dependent coupling coefficient 7. The 7^^ and 7e are not necessary equal but if we 
equate the absorption rates from one photon Fock states (i.e. |1)) we obtain that 7^^ = 7e. 

The one-count operator (Eq. (|T2l) ) and the probability (Eq. dHl)) and the no-count 
operator (Eq. ( |T3l) ) and the probability (Eq. ( |T5l) ) given by our model are equivalent with 
those given by the SD model (see App. IA2I and IA4I) . Thus our model gives an ab-initio 
derivation for the initially postulated SD model. 

Table [T] summarizes the expectation values of the number of photons after the one-count 
event, reported previously in references [g], [l^. For detailed derivations see Appendix lA 3[ 
Also the probabilities of observing the vacuum state after the one-count event are given. Note 
the difference between the models for the thermal and coherent ffelds. Our model agrees for 
all calculated expectation values with results given by the SD-form of superoperators while 
they disagree with the results given by the E model as seen in Table HI 

TABLE I: The expectation values of the number of photons after a one-count event are given in 
the upper part. In the lower part the probabihties of the vacuum state after one-count event are 
given for the thermal and the coherent fields (po and pi are the probabilities of the zero and one 
photon number states, respectively). If the field is initially in the Fock state |1) the probability 
Poit^) = 1 after the one-count event, otherwise po{t^) = 0. 



Initial state 


SD 


model 




E model 


Fock 

Thermal 

Coherent 


n(t+) 
n{t+) 
n{t+) 


=n{t) - 1 

=2n{t) 

=n{t) 


n(t+ 
n(t+ 
n(t+ 


)=n(t) - 1 

)=m 


Thermal 
Coherent 


Po(t+) 
Poit+) 


1 


Po{t^ 
Po{t^ 


)- 1 


-(l+n(t))^ 


l+n(t) 

\_ "{*) 
) e"(*)-l 



Table [m gives photon count rates for the SD and the E models. The results in Table ITTl give 
the probability w^^\t)dt of counting one photon at [t, t-l- dt] and the conditional probability 
w^^\t'^\t)dt of counting one photon at [t,t + dt] and the second photon at [t'^,t'^ + dt] 
[t^ = t + dt is infinitesimally greater than t and dt is the infinitesimal time that the one- 
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count process takes i.e. we are calculating the probability of absorbing the second photon 
immediately after the first photon). Furthermore, the second order coherence degree 

w«(t+|t) 



(18) 



is also given in Table [TTl For the derivation of the photon correlation function see Appendix 
lA 7[ Our results agree with the results of the SD model. 

TABLE II: The one-count rates {w^^^{t)), conditional one-count rates (w^^^f^lt)), and the second 
order coherence degrees {g^'^\t,t~^)) given by the SD and the E models. 



Initial state 



Fock 

Thermal 

Coherent 



lsdn{t) -Jsdinit) - I] 
-fsdn{t) 2-isdn{t) 
lsdn{t) 7sdn(t) 



7e 

n(t) 
7e l+n{t) 



7e(l 



n(t) 
7e l+n{t) 



gii(t)_l 



(2), 



n(t)-l 
n(t) 

2 
1 



1 
1 

e"(*)-(n(f)+l) 

gn{t)_|_g-n{t)_2 



The evolution of the expectation value of the number of photons in the SD and E models 
are respectively (see Appendix lA 6p 

fisdit) = n{0)e-^^'^' (19) 

ne(t) = n(0) +7e / iPoit') - l)dt'. (20) 
Jo 

Thus n{t) decays exponentially in the SD model irrespective of the field type while in the E 
model n{t) depends on the time integral of the vacuum state probability which is different 
for different field types. Again our results agree with the results of the SD model. 



V. DISCUSSION 

The E model of photon counting was originally motivated by arguing that the results 
given by the SD model are unphysical. Next we will show by selected examples that, on the 
contrary, it is the SD model that gives physical and also intuitionally right answers. 

In references Dodonov et al. gave the following main arguments in favor of the 

E model: (1) The expectation value of the number of photons may increase when operating 
with the SD one-count operator 



101 ]. see Table [B (2) The absorption rate of photons in the 



SD model is proportional to the number of photons and does not saturate for high n. 



1. Increase of n after one-count event 

In the SD model the expectation value of the number of photons after the observed 
absorption of one photon (i.e. after the one-count event) may be greater than before the 
absorption. For example after detecting of a photon with the SD one-count operator from 
the thermal field the expectation value of the number of photons doubles (see Table [T]). The 
assumption that this feature of the SD model is unphysical was the main justification for 
the E model. This behavior caused by the measurement back action to the cavity field 
is, however, in agreement with the photon bunching effect: When a photon arrives to a 
detector the expectation value of the number of photons raises and thus it is more probable 
to detect another photon. The growth of the expectation value of the number of photons 
is not unphysical and means that the states which had one or more photons become more 
probable after the absorption of the first photon. 

Furthermore, we have shown that in our model (as well as in the SD model) the expec- 
tation value of the number of the photons has exponential decay in time for all fields as it 
should. Thus, even if the expectation value of the number of photons may increase in the 
one-count event, on average the expectation value of the number of the photons decreases 
on every time interval during which no measurement induced projection of the field state 
occurs. 

2. Saturation of absorption rate 

The one-count rates for the SD and E models are wf}{t) = 7sdn(t) and wi^\t) = 7e(l — 
Po{t)), respectively. The fundamental difference between the models is that in the SD model 
the photocount rate is proportional to the expectation value of the number of photons while 
in the E model it is proportional to the probability that photons exists. Thus the absorption 
probability saturates in the E model. This may be a reasonable assumption for the very 
high intensities when the detector itself may saturate. However, the possible saturation 
of absorption probability is a internal property of the detector and should be an auxiliary 
model property, not a built in property of the quantum jump superoperators. In principle 
it is possible that at high intensities the detector could scatter some photons back to the 
cavity. This would lead to non-constant 7 which the present models (SD or E) can not adopt 
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as such. 

Let us consider two simple examples: (a) The field is in the state |^) = (|0) + |l))/v^. 
The one-count rates are w^J'J{t) = 7sd/2 and wi^\t) = 7e/2. (b) The field is the state |^) 
= (|0) + |100))/V2. The one-count rates are w^i^(t) = 1007,^/2 and w^e\t) = 7e/2. The 
behavior of the SD model seems more reasonable because it should be more probable to 
detect photons when the expectation value of the photons is larger. 



3. Second order coherence degree 

The photon bunching effect is related to the second order coherence: If < ^''•^''(O) < 1 
the light is antibunch ed, if g^'^\0) = 1 the light is nonbunched or random, and if g^'^\0) > 1 
the light is bunched Here the argument zero means that the correlation is measured 
with zero time delay. Examples of light fields obeying the above conditions are the Fock 
state, the coherent state and the thermal state, respectively. The calculated second order 
coherence degrees are given in Table [Tll 

For single-mode fields the second order coherence degrees are independent of the time 
delay between the measurements. The values are g^"^^ = 1 for the coherent field, g^'^^ = 2 for 
the thermal field and (yf^^-* = m/ {m — 1) for the Fock state \m) (m > 2) \v\. Comparing the 
results given in Table [III we note that the SD model reproduces the second order coherence 
degrees of the three standard example fields. 

The results given by our model (as well as the SD model) are equal with the experimental 
results given in Reference 8| for thermal and coherent fields. Note, however, that the 
measurement results in Reference js] are not given for light fields in a closed cavity but 
for continuous wave laser with a Gaussian frequency distribution for thermal light. We are 
comparing the results measured with zero time delay since the zero time delay cancels the 
multi-mode effect. Furthermore, we expect that the measurement results of the second order 
coherence degree in free space correspond to our theoretical results since g^'^^ depends only 
on the statistics of the light field, not the spatial distribution. 

n 

Dodonov et al. [6| also concluded that the SD model always predicts the photon bunching 
phenomenon for any initial field. We have shown (see Table [Tll) that this conclusion given 
in reference M is incorrect. 
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VI. CONCLUSIONS 



In this paper we have shown that by starting from the time- dependent Schrodinger equa- 
tion one can derive an exact difference equation for the density operator p/ of the combined 
detector-field system. Reducing this change of the density operator with respect to the de- 
tector variables gives a time integrable master equation for the field density operator. This 
gives the first principle expressions (Eqs. ffT^ and ffT^ ) for the one-count and no-count 
operators. Our results for the one-count and no-count operators agree with the SD model 
(see App. IA2l and IA4l for comparison) which therefore, in spite of being initially introduced 
ad hoc [4|, is actually an exact quantum mechanical result. 

APPENDIX A: DERIVATIONS 
1. Derivation of master equation 

The time-dependent Schrodinger equation is 

i^^l^.W) = (^o + ^;)I^.W), (Al) 

where = hwoo^a + fujjA\e) (e| is the time-independent field + detector Hamiltonian and ff| 
= hVt [a\e){g\ + a)\g){e'\) is the time-independent interaction Hamiltonian. The interaction 
representation is given by unitary transformation 

|vl>,(t)) = e'^o*/'^|vi>,(t)), (A2) 
where |\E'/(0)) = |\E's(0)). The Schrodinger equation gives 

= -i^o1^/(^)) + e'^o*/^(i^o^ + H'i)e-'^^'l^\^ i{t)) 

= Hn^iit)), (A3) 

where the last form holds, since [Hq, Hf\ = 0. Note that, since Hq and Hf commute, the 
interaction Hamiltonian in the interaction picture is same as in the Schrodinger picture. 
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Equation (1A3P gives following differential equation for the density operator of the coupled 
field-detector system in the interaction picture 



ih- 



dt 



The density operator can be written using the Taylor series as 

dp/(t) 



pi{dt) = p,(0) + 



dt 



t=0 



{dtf + .... 



(A4) 



(A5) 



2 dt2 

Since the interaction Hamiltonian is time independent, equation (1A4I) can be used to obtain 
the higher derivatives in the Taylor series 



dV/(t) 

dt2 



L 



Taking the terms up to second order in dt the Taylor expansion can be written as 

p{dt) = p,(o) - ^[i^;,p/(o)]dt - ^ [hi [Hipim] m\ 



(A6) 



(A7) 



(A8) 



where the double commutator is 
Pi{0)HfHf - 2Hfp{0)Hf giving 



Hf,[Hf,pr 



[HI H^jpriO)-pr{0)Hf] = ^|i^|p/(0) + 



p{dt) = pj{0)-^[HfJj{0)]dt+^{2H^jmHl-{H^jH%pj{0)} 



(A9) 



2. One-count probabilities and rates 



The one-count probability in the time interval [t,t + dt] is given by 7ATrace{yi^y4p/}dt. 
Thus 7^Trace{/i^Ap/} is the photon count rate. For the SD and E models we, respectively, 
obtain the photon count rates ( Trace{-} = Yl'^=o{^\ ' \^)) 



Wsd{t) = Isd npnit) = 'jsdnjt) 

71=0 

oo 

We{t) = TeX^PnW = 7e(l -Po(i)) 



(AlO) 
(All) 



n=l 



The one-count probabilities are obtained by multiplying the rates with dt. Thus 



one— count 



it) 



%dn{t)dt and 



one— count 



(t) =7e(l-po(t))dt. 
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3. One-count event 



If the one-count event is detected the density operator must change in accordance to the 
operation by the one-count operator and normahzation. The density operator after one- 

count is pj{t+) = ^^^^igf,)^t} = Tvact{lt%(t)} - and E models we, respectively, 

obtain 

^ CO 

^ ' n,n'=0 



n{t) 



J2 Pn+i,n'+i (t) V(n + l)(n' + l)|n) {n' \ (A12) 



n=0 



3-—- p^,At)\n-l){n - 1| = ^^-j- . (A13) 



n,n'=l 

Thus the new probabilities of the field states of n photons are 



■n -I- 1 



4. No-count probability 



The no-count operator is SrPf = ef^'^l^pfe^^'^l^^ where Ya — — i-^o — \ia^^A. and 
^0 = hujQa^a. The series expansion in the SD model is 

e{---^-h¥'-\m) = ' ^^r^'''''^""V) - e(---^-)-nm) (A16) 

n=0 

Thus the time evolution of the density operator in the no-count event is 

00 

m,m'=0 

The series expansion in the E model is 

cx 



^ ,-ia;o)"(ata)"r" ~ {-\^)\E^EfT^ 

giving 



2 

n=0 fe=0 



e^-|m) = e-'"""*" ^ ^'==° ^ ' ^ { ' ^' (A18) 

|0) |0), m = 



14 



and the evolution of the density operator 

oo 



m,m'=0 

Po,o 



|0)(0| + J2 kn^|0)(m|e+^-o"^--^- +p^,o|m)(0|e-^-«^'"--^^ 



m=l 



+ J2 Pm,m'e-'"^°^"'^^-^^\m){m'\. (A19) 



m,m'=l 



We are interested only on the diagonal elements to be able to calculate the evolution of the 
probability of the n photon state by taking the matrix elements {n\ ■ \n). Thus the SD model 
gives 

(n|5',p/(0)|n) = e-'^^->„(0) (A20) 

and the E model gives 

{n\SrPfmn)={ (A21) 
\ Po(0), n = 

The no-count probabilities at the time interval [t, t + r] are the sums over the probabilities 
{n\SrPf{t)\n) . The SD model and the E model, respectively, give 

oo 

P;/-'=°™*(t,t + r) = Ye-''^^^^pn{t) (A22) 

n=0 

pno-count^^^ t + t) = p,{t) + (1 - Po(t))e->^ (A23) 



5. No- count event 

If we know that the one-count event has not happened, the density operator must change 
in accordance to the operation by the no-count operator and normalization. Thus the 
probabilities of n photon states are obtained by normalizing equations flA20p and flA2ip 
giving for the SD and E models, respectively, 

Pn{t + r) = ' (A24) 



p^{t + T) = I foW+{i-PoW)e--- - (A25) 



n > 



po(t)+{l-po(t))e-Te^' 



-, n = 0. 
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6. Evolution of expectation value of photon number 

The density matrix evolves according to equation 

^ = -iuj [a^ap - pa^a) + (-fAApA^ - ^(iUp + pA^A)) . (A26) 

The probabihties of n photon states are given by the diagonal elements {n\ ■ \n). Thus we 
obtain for the SD model 



dt 

and for the E model 

dpn>i{t) 
dt 

dpo(t) 



{{n + lK+i(t) - pn{t)n) (A27) 



dt 

The expectation value of the number of photons is defined as 



7e(p„+i(t) -p„(t)) (A28) 
lePiit). (A29) 



n=0 

For the SD model the time derivation of is given by 

n=0 n=0 



lt) = J2npn{t). (A30) 



oo 



7. 



( X]*^^ + l)npn+i{t) - y^n^pn(t) I = 7sd I - l)p„(t) - ^n^Pn(t) 

vn=0 n=0 / \n=l n=l 

oo 



n=0 

This has a solution 

n(t) = n(0)e-^^'^*. (A31) 

For the E model we correspondingly obtain 

dnjt) ^ dpnjt) ^ 

= 2^ ^~;if~ =^^1^ [npn+i{t) - npn{t)] 

n=0 n=0 



oo oo \ / oo 



7e 



^npn+i{t) - ^np„(t) j = 7,1 ^(n - l)p„(t) - ^np„(t) 

vn=0 n=0 / \n=l n=l 

oo 



-7. 



n=l 



Time integration gives for the expectation value of the number of photons a solution 



n{t) = n(0) + 7e / (po(t') - l)dt'. (A32) 
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7. Second order coherence degree 



The second order coherence degree is ll|, [12 



g (r„ti,r„t„r„t„ri,t,) - r^, tOGa)(r„ t„ r„ t,) ' ^^^^^ 

where 

G'(^)(ri,ti,r2,t2,r2,t2,ri,^i) 
= Tr{p^E(-)(ri,ti)E(-)(r2,t2)EW(r2,t2)...EW(ri,ti)} 

with E^~\r,t) and E(+)(r,t) being the negative and positive frequency parts of the electric 
field operator. The two-fold delayed coincidence rate i.e. the counting rate per (unit time)^ 
is given by [3] 



w 



(^)(ri,ti,r2,t2,r2,t2,ri,ti) = s^G(^)(ri, ti, r2, ^2, r2, ^2, ri, i^i), (A34) 



where s is the sensitivity of the detector. We consider only the temporal correlation so we 
assume that all of the position vectors are equal and drop the spatial coordinate. We can now 
use well know formula of conditional probability: the probability that an event B occurs with 
the condition that A has happened is p(5 1 A) = p{BnA)/p{A). Th.usp{BnA) = p{B\A)p{A) 
giving w^'^\t'^ ,t){dtY = w^-^\t'^\t)dt w^^\t)dt, where we are considering correlation with 
infinitesimal time difference. Furthermore, we can write the second order coherence degree 
using the count rates 

where we, furthermore, assume that w^^\t~^) = w^^\t) due to the differential time difference. 
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